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00 \ We compare, for the static potential and at short distances, perturbation theory 

■ with the results of lattice simulations. We show that a renormalon-dominance picture 

. explains why in the literature sometimes agreement, and another disagreement, is 
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^ • perturbation theory agrees with lattice simulations. 
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1 Introduction 



The static potential is the object more accurately studied by (quenched) lattice simulations. 
This is due to its relevance in order to understand the dynamics of QCD. On the one hand, it 
is a necessary ingredient in a Schrodinger-like formulation of the Heavy Quarkonium. On the 
other hand, a linear growing behavior at long distance is signaled as a proof of confinement. 
Moreover, throughout the last years, lattice simulations [1, 2, 3] have improved their predic- 
tions at short distances allowing, for the first time ever, the comparison between perturbation 
theory and lattice simulations. Therefore, the static potential provides a unique place where 
to test lattice and/or perturbation theory (depending on the view of the reader), as well 
as an ideal place where to study the interplay between perturbative and non-perturbative 
physics. This is even more so since the accuracy of the perturbative prediction of the static 
potential has also improved significantly recently [4, 5, 6, 7]. 

Let us first review the status of the art nowadays. The prediction for the perturbative 
static potential at two loops [4] indicated the failure (non-convergence) of perturbation theory 
at amazingly short distances. This failure of perturbation theory is not solved by the inclusion 
of the leading logs at three loops computed in [5, 6] nor by performing a renormalization 
group (RG) improvement of the static potential at next-to-next-to-leading log (NNLL) [7].^ 

On the other hand, it was soon realized that the static potential suffered of renormalons 
[10] and that the leading one (of O(Aqcd)) cancelled with the leading renormalon of twice 
the pole mass [11]. Nevertheless, the prediction of perturbation theory for the slope of the 
static potential should not suffer, in principle, from this renormalon and could be compared 
with lattice simulations, which, indeed, only predict the potential up to a constant. This 
comparison was performed in Ref. [12], where they compared the RG improved predictions 
(without including ultrasoft log resummation) of perturbation theory up to next-to-next-to- 
leading order (NNLO) with lattice simulations. They indeed found that the discrepancies 
with lattice simulations and the lack of convergence of the perturbative series remained. They 
also found that the difference between perturbation theory and lattice could be parameterized 
by a linear potential in a certain range. Therefore, it seemed to support the claims of some 
groups [13, 14] of the possible existence of a linear potential at short distances. Such claims 
contradict the predictions of the operator product expansion (OPE), which state that the 
leading non-perturbative corrections are quadratic in distance at short distances^. 

On the other hand, driven by the analysis of Refs. [15], it has been argued [16] that, 
within a renormalon based picture, perturbation theory agrees with the phenomenological 
potentials aimed to describe Heavy Quarkonium. For his study, he used the static version 
of the IS mass and the upsilon expansion [17], which cancels the leading renormalon and 
relate the IS mass to the MS mass. Moreover, In Ref. [18] (see also [16, 19]), it has been 
shown that perturbation theory can indeed reproduce the slope of the static potential given 

^ There also exists a computation of the running of the Coulomb potential in vNRQCD [8] with the same 
precision that disagrees with the one obtained in pNRQCD [7]. At this respect, wc would like to report on 
a recent computation [9] of the 4-loop double log term of the Coulomb potential proportional to C^/3o that 
agrees with the pNRQCD result and disagrees with the vNRQCD one. 

^This contradiction is indeed so only if one believes that perturbation theory can be applied for the 
shortest distances available in lattice simulations ~ 1 — 4 GeV. We thank V.I. Zakharov for stressing this 
point to us. 
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by lattice simulations at short distances by using the force instead of the potential as the 
basic tool without the need to talk about renormalons. 

In this paper, we would like to try to clarify further the above issues and support the 
renormalon dominance picture by comparing the potential computed in lattice with the RS 
static potential (a renormalon free definition of the potential obtained in Ref. [20]). 

2 On-shell scheme 

The energy of two static sources in a singlet configuration reads (in the on-shell (OS) scheme) 

£;(r)=2mos+ lim ^H^u) ■ (1) 

1 — >oo 1 

In the situation Aqcd ^ ^l^i it can be computed order by order in a^iy) (i^ ~ 1/r) and in 
the multipole expansion (see [21, 7]) 

E{r) = 2mos + l^s,os(r, i^us) + (^^us(r, vvs) ■ (2) 
mos is the pole mass. The static potential reads {pts = ots{v)) 

oo 

V;,os = Vs,os{r, uvs) = E ^-^r'' (3) 

n=0 

where = 14(r, i/, i^us). The first three coefficients Vq, Vi and V2 are known [4]. The 
log-dependence on v of can also be obtained by using the i/-independence of T4,os- The 
log-dependcncc on i/us of V3 is also known [5] . Therefore, the only unknown piece of V3 is a 
z//z/us-independent constant. Its size has been estimated in Ref. [20] assuming renormalon 
dominance. For n/ = 0, it reads 

V^{r, 1/r, 1/r) = 1/r x (-76.1075) . (4) 

Leaving aside the leading-renormalon independent contributions to I3, this number suffers 
from two sources of errors (see Ref. [20] for details) . One is the error in the evaluation of the 
normalization factor of the renormalon, Nm-i which we will discuss in the next section. The 
other error is due to 0{l/n) corrections (n = 3 in this case), which appear to be negligible. 
The relative 0{l/n) corrections are -0.0264977, and the 0(l/n^) corrections are -0.00544012. 
Therefore, we expect the unknown 0{l/in?) effects to be very small. We will use Eq. (4) in 
the following for our estimates of V^,. 

It is quite remarkable that the numbers obtained in Ref. [20] for V3 are quite close to 
the ones obtained in Ref. [22] using Pade-approximants methods. This is even so if we do 
the Fourier transform back to momentum space. We display in Table 1 the numbers one 
obtains from Ref. [20] for the coefficient cq (the portion of the three-loop momentum space 
coefficient of the static potential independent of scale logarithms as defined in Ref. [22]) for 
different number of massless flavours. We can see that the difference with Chishtie and Elias 
results (see table 1 in Ref. [22]) is always of the order of 20-30 and roughly independent of 
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Table 1: Values of cq, as defined in Ref. [22], according to the results of Ref. [20]. 

the number of flavours. Therefore for ^ 0, where the absolute value of Cq increases, the 
bulk of Co agrees with Chishtie and Elias prediction. 

The static potential is i/-independent. By setting u = 1/r, we could effectively resum the 
In ur terms. The RG-improved expressions would read 

oo 

K,os = E Vn{r, 1/r, i/us)as(l/r)"+S (5) 

n=0 

and we would have expressions for n = 0, 1, 2, 3. In the above expression we have considered 
In rz/us ~ 1- Since In rz/us 3> 1 for some range of the parameters, one could be in the situation 
where one also has to resum ultrasoft logs as it has been done in Ref. [7]. Nevertheless, 
explicit calculations show that, at least for the range of parameters studied in this paper, 
higher order ultrasoft logs are subleading even if sizable. However, for definiteness, we will 
work with the rcsummcd expression (the physical picture does not change anyhow), which 
we will add to V^lr, 1/r, 1/r): 

Vs{r, 1/r, u^s) ^ V^{r, 1/r, 1/r) + ^al{r-') log f ^4^) ■ (6) 

In principle, it could be considered to be more consistent to add the resummed expression 
for the ultrasoft logs to V2. Nevertheless, we will add them to V3 since it is the counting 
consistent when the leading single log gives the bulk of the correction (moreover, we will 
see that, even after the leading-renormalon cancellation is achieved, its effect is smalP if 
compared with the typical estimate of the terms of V;:;). 

The leading O(r^) non- vanishing contribution to SE\js{r,i'us) reads (we actually write 
the Euclidean expression) 

3iVc Jo 

1 as 

At the order of interest Vgos = ^ and Vg os is approximated to its leading order value: 

2Nc r 

as 

—Cf — • This contribution cancels the leading z/us scale dependence of Vg os- Foi' the scales 
r 

we will study in this paper (those accessible by lattice simulations), we will consider a^/r to 
be a non-perturbative scale. Therefore, Eq. (7) can not be computed perturbatively and it 
will not be considered in this and the following section, where we aim to a pure perturbative 
prediction (it will be considered, however, in sec. 4.1). Thus, our results will depend on i/^s, 
for which we will take Uus = 2.5 r^^. Wc will set n/ = and work in lattice units: r^^ (with 
rg^ ~ 400 MeV according to lattice simulations), — 0.602(48) r(7^ [23], since our aim 
will be to compare with lattice simulations. 
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Figure 1: Plot o/roVbs(r) at tree (dashed line), one-loop (dash-dotted line), two-loops (dotted 
line) and three loops (estimate) plus the leading single ultrasoft log (solid line). For the scale 
ofas{u), we set v = 1/0.15399 

''"O ■ ^us — 2.5 Tg . 

Now, we want to check the convergence of the series within perturbation theory. We 
first consider the perturbative static potential without trying to resum Inru logs, ie. for 
a fixed scale u for as(^) (Eq. (3)). We show our results in Fig. 1, where we have chosen 
u = 1/0.15399 rg'^. We see that it does not converge to a value. We next try with Inru 
resummation (Eq. (5)), since these logs could be large. We show our results in Fig. 2. We 
find the same problem and moreover the slope appears to be different in both cases. The 
situation appears to be puzzhng. 

Let us now compare with lattice simulations. These predict the energy of an static 
quark- ant iquark system up to an r-independent constant^: 



Therefore, the slope of the potential is a physical prediction of the lattice static potential 
and does not suffer from the leading renormalon. The constant is usually fixed in a way such 
that 



except in Ref. [2] where -Eiatt.l'^c) = is used (for the definition of Tc see [2]). 

In order to compare the perturbative predictions with lattice simulations, we add a 
constant to the perturbative potential to fix both potentials to be equal at some scale r': 



■^However, it remains to be seen whether, once V3 is exactly computed and the O(r^) renormalon were 
also subtracted, the ultrasoft logs would be a subleading effect compared with finite pieces. 

''Therefore, due to the unknown constant it can not be related with the underlying theory (QCD) and, 
for instance, obtain the mass. 



^iatt.(r) = fciatt. + ^lim^ ^ ln{Wn). 
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-Eiatt.('^o) = , 



(9) 




(10) 




Figure 2: Plot ofroVos{r) at tree (dashed line), one-loop (dash-dotted line), two-loops (dotted 
line) and three loops (estimate) plus the RG expression for the ultrasoft logs (solid line). For 
the scale of as{i^), we set u = 1/r. u^s = 2.5rQ^. 

The value of this constant varies according to the order the perturbative computation has 
been done, ie. 

fcpert.(r') = fco(r') + fcias(r') + k2al{r') + ■■■ . (11) 

We will consider the lattice data of Refs. [1, 2]. They both appear to be perfectly 
compatible with each other for finite lattice spacings. In any case, the lattice data of Ref. 
[2], being more recent, appears to have smaller errors with a very small spread of the points 
around the fitted curve. Here, in order to get the highest possible precision, we will use the 
lattice data of Ref. [2] for which the continuum limit has been reached. For comparison with 
the lattice points at finite lattice spacing of Ref. [1] and for a somewhat a larger range, we 
refer to the 1st version of this paper sent to the web. The physical outcome is, in any case, 
the same. 

The comparison with lattice is performed in Fig. 3 for the perturbative result without 
ln(ri/) resummation (z/ constant) and in Fig. 4 for the perturbative result with ln(rz/) 
resummation {u = 1/r). As explained above, we have added a constant to the perturbative 
potential (different at each order) as to make it equal to the lattice potential at r = 0.15399 tq. 
We see a completely different behavior in both cases. Whereas for the latter, the agreement 
goes worse as we increase the order of our calculation, for the former the agreement with 
lattice improves^. Indeed, it is quite remarkable that lattice simulations are precise enough 
to see the log behavior of predicted by asymptotic freedom. 

In the next section we will explain the above behavior (why sometimes perturbation 
theory can describe lattice data and some others not) within an controlled framework where 
the renormalon is taken into account. 

^Usually, the comparison between lattice and perturbation theory was performed using the expressions 
with V = 1/r, whereas the analysis with v constant was largely unnoticed. 
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Figure 3: Plot o/ro(Vbs(^) — Vos(^') + -E'zatt. (^') ) versus r at tree (dashed line), one-loop 
(dash-dotted line), two-loops (dotted line) and three loops (estimate) plus the leading single 
ultrasoft log (solid line) compared with the lattice simulations [2] EiattX''^)- ^^r the scale of 
as{i^), we set v = 1/0.15399 Tq"^. = 2.5rQ^ and r' = 0.15399ro. 



3 RS scheme 

The RS scheme was defined in Ref. [20] aiming to ehminate the renormalons of the matching 
coefficients appearing in Heavy Quarkonium calculations. In particular, the leading renor- 
malon of the heavy quark mass and the static potential. We refer to Ref. [20] for details. 
Here, we just write the relevant formulas needed for our analysis. 

Analogously to the OS scheme, the energy of two static sources in a singlet configuration 
reads 

Esir) = 2mnsiiyf) + f hm ^ \n{Wu) - 25mRs(^//)) , (12) 

where^ 

n=l n=l \^^/ k=0 iU + O-Kj 

and m^g = mos — <5^rs- If the beta function were known to infinity order in perturbation 
theory, it would be possible to obtain all the coefficients b and Ck [24]. In practice, only b 
and co,i,2 are known (see [24, 20, 25]). For Nm only an approximate calculation is possible 
using some ideas first developed in Ref. [26]. This computation has been done in Ref. [20]. 

^Actually, we are going to use in this paper the RS' scheme defined in Ref. [20] instead of the RS scheme, 
since we beheve it has a more physical interpretation. Nevertheless, we will denote it in this paper as RS 
scheme in order to simplify the notation. In any case, the physical picture does not change. The results 
would also converge towards the result predicted by lattice simulations but with a different pattern. 
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Figure 4: Plot of rQ{Vos{r) — Vos{r') + EiattX^')) versus r at tree (dashed line), one-loop 
(dash-dotted line), two-loops (dotted line) and three loops (estimate) plus the RG expression 
for the ultrasoft logs (solid line) compared with the lattice simulations [2] Eiatt.i''^)- ^^r the 
scale of as{h'), we set v = 1/r. u^s = '^^-^Tq^ and r' = 0.15399 tq. 



The result reads (for Uf = 0) 



Nm = 0.424413 + 0.174732 + 0.0228905 = 0.622036 , 



(14) 



where each term corresponds to a different power in u (in the Borel plane) of the calculation 
(see Ref. [20] for details). We see a nice convergence. This number will be the one we will 
use in the following. This number should be equal to —2Nv, where A'^v' is the normalization 
factor of the renormalon of the static potential. Ny was also approximately computed in 
Ref. [20]. For ri/ = 0, it reads 



Nv = -1.33333 + 0.499433 - 0.338437 = -1.17234 . 



(15) 



In this case the convergence is not so good but we have an alternating series. In any case, 
we see that both values appear to be quite close: 



2Nm + N^ 



0.059. 



(16) 



2N^ - Nv 

We take this as an approximate indication of the error in the evaluation of Nm- 

In the situation Aqcd ^/f, Es{r) can be factorized in the following way (with the 
accuracy we are working, the expression for 6Es,us will be equal in the OS and RS scheme): 

Es{r) = 2mRs(z//) + Vs,Rs{r; Vus\ Vf) + 5Es,vs{r, Vus) ■ (17) 
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V^3s(^) = V^,Rs('"; J^us', ^f) was defined in Ref. [20] (note that we expand Vs and 5m^s with 
the same ccg): 

oo oo 

K,Rs(r) = K + 2(5mRs - E + 2<^^RS,nX = E yiln<^'- (18) 

n=l ra=0 

Due to the leading renormalon, the behavior of the perturbative expansion at large orders 
of Vn reads [24, 20] 

Therefore, by adding 5mRs to Vs we expect to cancel the renormalon of the static potential. 
Indeed, what we want is to obtain the renormalon cancellation order by order in perturbation 
theory. This means to kill the 0{n\) behavior of the coefficient multiplying the 0(«") term 
of Vs with the 0{n\) behavior that appears in 5m-Rs such that Vrs enjoys nicer convergence 
properties. In principle, this cancellation should hold for an arbitrary u and, in particular, 
for i/ = 1/r. This, indeed, is the key observation of this paper and it will become relevant 
later on in order to explain the results in the previous section as well as in the literature. 

In principle, any i/f should cancel the renormalon of the static potential (we have the 
constraint, though, that Uf should be large enough to fulfill a;s(z/j) <^ 1). If we re-expand 
(5mRs(z/j) in terms of as(z/), the renormalon is still kept in 6m^s{'^f) but new terms are 
generated, which, nevertheless, should not belong to the renormalon. In other words 

Smjis{i^f) = SniRsii^) + F{u; Uf), (20) 

where the Borel transform of F{p; vj) should be analytic aXu — 1/2 in the Borel plane. This 
can be easily seen in the large (5q limit. In more formal terms what we have is 

^Ms ~ constant I m[5m-Rs,{i^f)] , ^ms ~ constant ImlSm-Rsi^)] 

with the same constant by definition (indeed this is what defines the renormalon). Therefore, 
the conclusion is that 

Im[F{iy;iyf)] = 

to the order of interest, concluding that F docs not have the closest singularity in the 
Borel plane. Note that in order to have the structure of Eq. (20) is necessary to have the 
complete renormalon contribution, ie. to all orders in perturbation theory. If we have the 
renormalon to some given order in perturbation theory, a change of the scale will produce 
large (rcnormalon-relatcd) terms that would not cancel since we do not have the complete 
sum. Indeed, this is the explanation why in fixed order perturbation theory one has to 
expand on ag at the very same scale for the renormalon and the OS result in order to achieve 
the renormalon cancellation order by order in ctg. 

On the other hand, Uf can not be arbitrarily large. Otherwise the power counting is 
broken. This restricts Uf < 1/r in order to ensure the counting Vrs ~ 0{as/r). In our case, 
we will set z/j = 2.5 Tq"^ ~ 1 GcV. In any case, when studying the slope of the potential, the 
result should be independent on Uf, as 5mRs(t//) is just an r-independent constant. At this 
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respect, we should say that this is indeed so to a large extent. If we take uj = 5rQ ~ 2 GeV 
agreement with lattice is also obtained^. One could also think of setting i/j = 1/r. This 
would mean to subtract an r-dependent constant to the static potential. This is not what 
it is done in lattice simulations since they are arbitrary up to a r-independent constant. 
Therefore, this would not agree with lattice simulations. 
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Figure 5: Plot o/roVRs(r) at tree (dashed line), one-loop (dash-dotted line), two-loops (dotted 
line) and three loops (estimate) plus the RG expression for the ultrasoft logs (solid line). For 
the scale of a^iy), we set v = 1/r. u^s = 2.5 r^^^ and z/j = 2.5 Tq^. 

In conclusion, by using Vrs, "we expect the bad perturbative behavior due to the renor- 
malon to disappear. Let us see that this is indeed so. Working analogously to the OS 
scheme, we consider the RS potential at different orders in perturbation theory with ln(rz/) 
resummation or not. We display our results in Figs. 5 and 6. In the first case, the ln(rz/) 
resummation is achieved by setting u = 1/r and the ultrasoft log resummation is also in- 
cluded. In the second case, we set z/ = 1/0.15399 and only the leading single ultrasoft log 
is included. We see that now we do not have the gap between different orders in the pertur- 
bative expansion (or it is dramatically reduced). Moreover, we do not have a (dramatically) 
different behavior of the slope of the potential; In(ri^) resummation or not converges to the 
same value. 

We can now compare with lattice simulations: Figs. 3 and 7. We work analogously to the 
OS scheme. Now, the constant we have to add is approximately ttg-independent, reflecting 
that we have accurately achieved the renormalon cancellation. We see that our results are 
convergent to the same potential, which, as we can see, corresponds to the lattice potential. 

^One could say that since Sm^s starts at 0{al) a larger value of could be used and yet not to break 
the power counting. Nevertheless, we refrain from using Vf ~ ^Tq^ since we do not understand the meaning 
of the results with values of Vf larger than 1/r. In particular, this would have a difficult explanation from 
an effective theory point of view, where Vf is understood as an cutoff of the effective theory, which has to 
be smaller than 1/r. 
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Figure 6: Plot o/roVasl?') at tree (dashed line), one-loop (dash-dotted line), two-loops (dotted 
line) and three loops (estimate) plus the leading single ultrasoft log (solid line). For the scale 
ofasii'), we set u = 1/0.15399 r^^ . u^s = 2.5rQ^ and uj = 2.5 r^^ . 

The case without ln(rz/) resummation is exactly equal to the OS case, ie. to Fig. 3, 
since it is just equivalent to add a r-independent constant at each order in Og. This explains, 
within the renormalon dominance picture, the success of that specific OS scheme calculation; 
it is just equivalent to an specific case of the renormalon based calculation: u = constant. 
Therefore, perturbation theory without log resummation can explain lattice data up to 
r ~ 0.5 ro by subtracting a r-independent constant at each order in perturbation theory. 
This fact is understandable if the factorization scale chosen is not very far of 1/r so that 
the RG does not play a decisive role plus the fact that subtracting a constant kills the 
renormalon. However, now, we are not restricted to take u = constant but we can take 
z/ = 1/r, which allows us to resum the ln(rz/) terms and yet cancel the renormalon as we can 
nicely see in Fig. 7. 

Within the renormalon dominance picture, we can also explain why perturbation theory 
with ln(ri/) resummation could not explain lattice data (see Fig. 4) even if subtracting a r- 
independent constant at each order in perturbation theory. In particular, this would explain 
the disagreement found in Ref. [12]. The explanation of this fact is that ln(rz/) resummation 
is equivalent to set as{r) as the expansion parameter. Therefore, the expansion parameter is 
a function of r and, consequently, the amount to be subtracted at each order in as{r) in order 
to cancel the renormalon is also r-dependent. It follows that the (perturbative) renormalon 
cancellation can not be achieved by only subtracting a r-independent constant at each order 
in as{r) as it was done in Fig. 4. 

We would like now to understand the agreement between lattice [18] (or potential models 
[16]) and perturbation theory if the static potential is reconstructed from the force in the 
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Figure 7: Plot o/ ro(VRs(r) — V^si^') + Eiatt.ij"')) versus r at tree (dashed line), one-loop 
(dash-dotted line), two-loops (dotted line) and three loops (estimate) plus the RG expression 
for the ultrasoft logs (solid line) compared with the lattice simulations [2] Eiatt.ij')- ^^r the 
scale of as{h'), we set v = 1/r. uj = Uus = 2.5 and r' = 0.15399 tq. 



following way: 



V{r) 



dr'F{r') + V{ro) 



ro 



ro dr 



(21) 



As far as we end up in the potential, we should have the same problems (and solutions) than 
in the previous sections. The point is how Eq. (21) is implemented in practice. The force 
also has an expansion in as- 



F{r) 



;i/r) 



+ 



■}■ 



(22) 



In practice, the force is introduced in Eq. (21) order by order in a^- If the expansion 
parameter is as(i^), the potential is reconstructed order by order in asiv) and we are exactly 
in one of the situations considered in this and the previous sections (Fig. 3). If the expansion 
parameter in the force is as(l/r), as it was made in Refs. [16, 18], we are in a new situation. 
The introduction of a;s(l/r) in the integral in Eq. (21) produces and infinite series in Q;s(l/r) 
and as{l/rQ). At this respect, it is not clear what the expansion parameter is and the 
systematics of this procedure. In any case, here, we are just interested to see whether the 
renormalon cancellation is achieved for these calculations. At this respect, the point is that 
the expansion parameter in the potential is not as{r). From the computational point of view, 
all the powers in Q;s(l/r) and as(l/ro) are considered to be of the same order. If we consider 
the leading order, one loop, approximation we have aqq = ag and l3{aqq) = — /5oas/(27r) and 
we obtain 



dr'F{r') 



'ro 




Poas{l/r) 

27T 



En! 

n=0 



271 



(23) 
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Both terms have the renormalon (in the large /3o approximation) so that they cancel each 
other. We expect that higher orders computations will reconstruct the QCD renormalon 
but yet the renormalon cancellation should remain at any order since the same (partial) 
piece of the renormalon appears for f dr' F^~^°°^'^ {r') than for dr' F'^~''"°^^ {r') . This ex- 
plains within a renormalon- dominance approach the convergence of computations using this 
framework even if no reference to renormalon was made for these calculations. 

It is worth stressing that within this renormalon-based scheme calculation, we can now 
account for the logs correctly by setting u = 1/r (see Fig. 7). On the other hand, we 
introduce some errors when fixing u = 1/r, since we are only able to perform an approximate 
computation of the renormalon, which do not exist in the calculation of the slope with 
u = constant. The two things compete with each other in the final accuracy of the result 
(compare Figs. 3 and 7). Indeed, once the renormalon cancellation has been achieved, the 
resummation of logs does not appear to be very important. At least, if we do not choose our 
scale very far from the typical scales at study. If we go to scales smaller than 2.5 , the 
calculation with u = 1/r appears to give better results (compatible with lattice, although 
with larger errors, up to amazingly long distances (~ 0.8 ro)) than the calculation with fixed 
u (see Fig. 8). This may signal that log resummation becomes important for these scales 
(note that we have chosen u = l/0.15399r(^^ as the starting point of the log evolution, for a 
smaller value of u an improvement may be obtained). It should be studied further whether 
perturbation theory can indeed describe lattice data in the regime 0.5 — 0.8 ro, since this 
would support the claims of Ref. [16] in this direction. 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 



a) r/ro b) r/ro 

Figure 8: Figs, a) and h) correspond to Figs. 3 and 7 respectively up to the value r = Llr^. 

The main result of this section is to show that all the different results found in the previous 
and this section as well as in the literature about the convergence (or not convergence) of 
the perturbative series can be explained by the presence of the renormalon. We have also 
shown that, with a renormalon-based scheme, agreement with lattice is obtained. The point 
is that we now know the constant to be subtracted and its dependence in a^, therefore, we 
can set as at the same scale and obtain the cancellation for any r. In this way, we obtain 
the log resummation and the renormalon cancellation. One potential problem is that other 
logs (ln(rz//)) appear (in the physical, finite mass, case this problem does not really show up 
because one usually sets z// < 1/r being r the typical size of the system). 



12 



4 Bounds on short-distance non-perturbative poten- 
tials 



We now study possible non-perturbative effects in the static potential. The first thing to 
notice is that any non-perturbative effects should be small and compatible with zero since 
perturbation theory is able to explain lattice data within errors. We can make this statement 
more quantitative by using the lattice data obtained in Ref. [2] where the continuum limit 
has been reached. For these lattice points the systematic and statistic errors are very small 
(smaller than the size of the points). Therefore, the main sources of uncertainty of our 
(perturbative) evaluation come from the uncertainty in the value of (±0.48 Tq^) obtained 
from the lattice [23] and from the uncertainty in higher orders in perturbation theory. We 
show our results in Fig. 9.^ The inner band reflects the uncertainty in whereas the 
outer band is meant to estimate the uncertainty due to higher orders in perturbation theory 
by allowing cq to change by ±146. This is half its value according to the estimate in Table 
1. This may seem a conservative variation if we take into account that the difference with 
the estimate obtained using Pade-approximats [22] is Scq ~ 20. 




0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 



Figure 9: Plot of roiVRsi^) — ^s('^') ± Eiatt.ij')) versus r at three loops (estimate) plus the 
leading single ultrasoft log (dashed line) compared with the lattice simulations [2] Eiatt.ij)- 
For the scale of otsiv), we set v = 1/0.15399 Tq ^. u^s = 2.5 Tq and r' = 0. 15399 tq. The 
inner and outer band are meant to estimate the errors in and cq. For further details see 
the main text. 



^We have performed the estimation of the errors with the results of Fig. 3 (ly — constant) , since they do 
not introduce errors due to the evaluation of the renormalon. Nevertheless, a similar conclusion had been 
achieved if the analysis had been done with the results of Fig. 7 {v = 1 jr) . 
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4.1 Bounds on a short-distance quadratic potential 

The non-perturbative effects slioufd be encoded witliin Eq. (7). In general, tliis equation is 
a convolution of two scales: Aqcd and ag/r. The explicit functionality is unknown since it 
involves the knowledge of the details of the non-perturbative dynamics. In the formal limit 
1/r ^ Aqcd ^ ots/f, it is possible to know the scaling of the potential in Aqcd and r based 
on dimensional arguments [27]: 

Tt? r°° J- 
5Eus(r, i^us) ^ I dt{gE'^{t)<l>{t, 0):^^ gE\0)) (uus) ■ (24) 

This behavior for the non-perturbative potential is expected on the basis of the OPE 
in that limit. 

The fact that we are able to put bounds on our pcrturbative evaluation of the potential 
and that the errors of the lattice simulations arc very small allows us to put bounds on the 
non-perturbative effects. This is of utmost importance since these non-perturbative effects 
produce the largest errors in the determination of the bottom mass using the T(15') mass. 
Unfortunately, even though lattice simulations of the static potential are very precise, we can 
not obtain strong constraints on the non-perturbative effects. The main source of error is 
due to the uncertainty in Aj^. Assuming a non-perturbative potential V ~ Ar^, we obtain 

|^|<0.6ro-^ (25) 

from the variation of Aj^. If a potential 0.6 r^'^r'^ were introduced as a correction in the 
binding energy of the T{1S), it may induce a correction of ~ 100 MeV in the mass of the 
T(15') (however, let us also stress that lattice simulations are also compatible with A — 0). 
Nevertheless, if lattice determinations of A^^ could be improved, they would provide one 
of the most promising ways to control the errors in the determination of the bottom MS 
mass using the T(IS') mass. Let us also note that one can turn the problem around. If one 
assumes the non-perturbative effects to be small one can improve the lattice determination 
of Aj^ from the static potential data. 

4.2 Bounds on a short-distance hnear potential 

The usual confining potential, SV = ar, goes with an slope a = 0.21GeV^. In lattice units 
we take: a — 1.35 Tq^. Can we discriminate such potential at short distances with the 
available lattice data? The answer is yes as we can see from Fig. 10.^ The introduction of a 
linear potential at short distances with such slope is not consistent with lattice simulations. 
This is even so after the errors considered in Fig. 9 have been included. This should not 
come as a surprise since this linear potential appears as an effect of long distance (not short 
distance) and it is not expected to appear from Eq. (7). Therefore, it follows that the use of 
the Cornell potential (with the perturbative static potential emanated from QCD instead of 
a pure 1/r potential times a fitted constant) as a phenomenological fit of the static potential 
introduces systematic errors if the typical inverse Bohr radius scale of the heavy quarkonium 

^We have used the results from Fig. 3, since they do not introduce errors due to the evaluation of the 
renormalon. Nevertheless, the same conclusion had been achieved with the results from Fig. 7. 
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system to study lye in the short distance regime as it is the case, for instance, for the T(IS'). 



0.5 - 



- 



-0.5 - 



-1.5 - 

0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 



r/ro 



Figure 10: Plot o/ro(Vas('") — l^s(^^') + -^;att-('^')+'^('^~''^')) versus r at three loops (estimate) 
plus the leading single ultrasoft log compared with the lattice simulations [2] Eiatt.i''^)- ^or 
the scale of as{h'), we set v = 1/0.15399 Tq'"'^. z/„s = 2.5rQ^ and r' = 0.15399 tq. The dashed 
band is meant to estimate the combined error due and cq (see Fig. 9). 



On the other hand, recently, there have been claims about the possible existence of a 
linear potential at short distances [13, 14]. Expected it to be of different physical origin 
than the long distance linear potential, it may have a different slope than the (long distance) 
static potential discussed above. It would be very important to discriminate its existence, 
since such behavior at short distances is at odds with the OPE. This is indeed possible, 
since the short-distance linear potential expected in Ref. [13, 14] have an slope of the order 
of magnitude of the long-distance confining potential that we have already ruled out above. 
Therefore, we can conclude that no linear potential exists at short distance (with the present 
estimates for its slope). 



5 Conclusions 

In this paper, we have shown (Figs. 3 and 7) that perturbation theory can explain lattice 
simulations of the static potential up to scales of order ~ 2r^^ once the renormalon has 
been taken into account without the need of non-perturbative effects. Quite remarkable, 
lattice simulations are precise enough to see the scaling log behavior predicted by asymptotic 
freedom and not just the pure 1/r behavior of the Coulomb potential. 

We have shown that the different results one may find in the literature as well as in 
sections 2 and 3 about the convergence (or not) of the perturbative series to the lattice 
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result can be explained by the presence of the renormalon. In particular, we can explain 
why a calculation (Fig. 3) within perturbation theory can agree with lattice data, if the 
perturbative expansion is made with q;s(z/) with v constant, without an explicit reference to 
the renormalon, but not if the expansion is made with a^iXjr) (Fig. 4). This explanation is 
given in a renormalon-based approach in a unified way. Moreover, our method allows to set 
V —Ijr and, therefore, resum the \vl{tv) logs besides of canceUing the renormalon (Fig. 7). 

We have tried to find bounds on possible non-pcrturbative potentials at short distances. 
As a matter of principle, large non-perturbativc effects arc ruled out with the above analysis. 
In particular, they are compatible with zero (Fig. 9). We have shown (Fig. 10) that 
lattice data disfavours a linear confining potential at short distances with the usual slope 
(cr — 1.35 r(7^) assigned to the confining potential at long distances. For a possible non- 
perturbative potential, the main source of error comes from the one associated to the 
determination of in the lattice (the error due to higher order in perturbation theory 
appears to be less important). Therefore, we can not put very precise constraints on it in 
order to reduce the non-perturbative error on the determinations of the bottom MS mass 
from the T(15') mass. Nevertheless, lattice simulations remain as one of the most promising 
possibilities to constraint the size of the non-perturbative effects at short distances. 
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